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Abstract 

Fast spinning magnetars are discussed as strong sources of high energy neutrinos. 
Pulsars may be born with a short rotation period of mihiseconds with the mag- 
netic field amplified through dynamo processes up to ~ 10^^ — 10^^ G. As such 
millisecond magnetars (MSMs) have an enormous spin-down power ~ lO^'^ergs"^, 
they can be potentially a strong, extragalactic high-energy neutrino source. Specif- 
ically, acceleration of ions and subsequent photomeson production within the MSM 
magnetosphere are considered. As in normal pulsars, particle acceleration leads 
to electron-positron pair cascades that constrains the acceleration efficiency. The 
limit on the neutrino power as a fraction of the spin-down power is calculated. It 
is shown that neutrinos produced in the inner magnetosphere have characteristic 
energy about a few xlOO GeV due to the constraint of cooling of charged pions 
through inverse Compton scattering. TeV neutrinos may be produced in the outer 
magnetosphere where ions can be accelerated to much higher energies and the pion 
cooling is less severe than in the inner magnetosphere. High energy neutrinos can 
also be produced from interactions between ultra-high energy protons accelerated 
in the magnetosphere and a diffuse thermal radiation from the ejecta or from the 
interaction region between the MSM wind and remnant shell. The detectability of 
neutrinos in the early spin-down phase by the current available and planned neu- 
trino detectors is discussed. 
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1 Introduction 



Young pulsars are generally considered as an important high-energy neutrino 
source (e.g. Berezinsky & Prilutsky 1978; Sato 1978; Protheroe, Bednarek, & 
Luo 1998; Nagataki 2004). Particles can be accelerated either in the magne- 
tosphere by a rotation-induced electric field (e.g. Harding & Muslimov 1998; 
Hibschman & Arons 2001, and references therein), or in the pulsar wind, e.g. by 
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large amplitude waves (e.g. Gunn & Ostriker 1971; Sato 1978; Asseo, Kennel, 
& Pellat 1978; Usov 1994; Melatos & Melrose 1996) or magnetic field recon- 
nection (e.g. Lyubarsky & Kirk 2001), and interact with thermal photons from 
the star's surface or protons in the remnant to produce high energy neutrinos. 
There is strong observational evidence for pulsars with an extremely strong 
magnetic field exceeding the quantum electrodynamics (QED) critical field 

5 > 5c ~ 4.4 X 10^''^ G, called magnetar (e.g. Woods & Thompson 2004). Soft 
gamma-ray repeaters (SGR) and anomalous X-ray pulsars (AXP) are believed 
to be slowly rotating magnetars (Thompson & Duncan 1995). The possible 
existence of fast spinning magnetars has also been proposed (e.g. Usov 1994; 
Blackman & Yi 1998; Wheeler et al 2000; Rees & Meszaros 2000; Gaensler 
et al. 2004). Such MSM may be formed from an accretion- induced collapse, 
in which the magnetic field is amplified exponentially to 10^^ — 10^^ G by 
dynamo amphfication (Dar et al. 1992; Duncan & Thompson 1992; Thomp- 
son 1994). If such fast spinning magnetars exist they are by far the most 
energetic pulsars with a typical spin-down power 10^^ — 10^° erg s~^ for a mag- 
netic field Bq = 10^^ G and a rotation period P = 1ms, which is about 10^^ 
times more powerful than the Crab pulsar. MSMs spin down rapidly, produc- 
ing powerful transient gamma-ray emission within a typical spin-down time 
~ (1/2)7Qo/Le ^ 10^ s, where Qq and I ^ lO^^gcm^ are the initial angular 
velocity and inertial moment of the pulsar, and Le is the spin-down luminosity 
due to the magnetic dipole radiation. As in normal pulsars, rotation-driven ac- 
celeration leads to gamma-ray emission and in particular, acceleration of ions 
may lead to photopion production through either proton-photon or proton- 
proton interactions, producing neutrinos. Thus, MSM can be a potentially 
observable source of high-energy neutrinos during its initial spin-down phase. 

Neutrinos originating from magnetospheric acceleration has been considered 
by several authors for normal young pulsars (e.g. Sato 1978; Protheroe, Bednarek, 

6 Luo 1998; Bednarek & Protheroe 1997; Bednarek & Protheroe 2002). Neu- 
trinos from slowly rotating magnetars were discussed recently by Zhang et al. 
(2003). However, as such slow rotators have a much lower spin-down power 
than the Crab pulsar (by a factor of 10^), the corresponding luminosity of 
rotation-powered neutrino emission is marginal even for a km scale neutrino 
detector. So far, in most of the magnetospheric-origin models for neutrino 
emission, the constraints by pair production and radiation reaction were not 
considered. Magnetospheric acceleration arises from a rotation-induced electric 
field in the open field line region due to that the charge density of outfiowing 
charged particles deviates from the corotation charge density, referred to as 
the Goldreich-Julian (GJ) density. Pair cascades tend to screen out the ac- 
celerating electric field and send a backflow of opposite-charged particles that 
can reverse the sign of the electric field. Therefore, pair cascades can strongly 
limit the fraction of the spin-down power going into particle acceleration and 
constrain the neutrino luminosity. There are extensive discussions on neu- 
trino emission resulting from particle acceleration in the pulsar wind or at the 
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wind termination shock (e.g. Berezinsky & Prilutsky 1978; Sato 1978; Beall 
& Bednarek 2002; Granot & Guetta 2003; Nagataki 2004; Luo 2005). In these 
models, it is generally hypothesized that protons (or ions) can be accelerated 
efficiently to ultra-high energy. Although acceleration in the pulsar wind can 
be efficient in the sense that most of the Poynting flux that is thought to be 
dominant near the pulsar is converted to particle kinetic energy, as supported 
by the observational evidence from the Crab nebula (Gaensler et al. 2002), 
the specific acceleration mechanism is poorly understood (e.g. Kirk 2005) and 
there is no rehable estimate of the maximum energy of accelerated protons. 

In this paper, we explore the magnetospheric origin of neutrinos from MSMs in 
which, compared to normal young pulsars, acceleration processes are strongly 
modified by both the supercritical magnetic field and rapid rotation. There 
are basically two classes of model for particle acceleration in the pulsar magne- 
tosphere: the polar gap model, in which acceleration is assumed to occur near 
the polar cap (PC), and the outer gap model, in which the acceleration region 
is located in the outer magnetosphere. There is also a variant between the 
two, called the slot gap model (e.g. Arons 1983; Muslimov & Harding 2004). 
We consider these two cases separately and comment briefly on the slot gap 
model in Sec. 6. For acceleration near the PC we extend the polar gap model 
for normal pulsars (e.g. Hibschman & Arons 2001; Harding & Muslimov 1998) 
to MSMs. Since the steady gap model may not be realistic as pair cascades 
are likely nonstationary and time-dependent, we emphasize the energetics of 
polar gap acceleration rather than a specific gap model. The energetics, de- 
scribed by the acceleration efficiency, which is defined here as the ratio of 
pair-production limited potential to the maximum potential (across the PC), 
determining how much of the spin-down power goes into particle acceleration, 
can be determined from the relevant free path for pair production. A pair 
may be produced through interaction of a Lorcntz boostered thermal photon 
with the Coulomb field of an ultrarelativistic ions. The more efficient pair 
production process is single photon decay in a strong magnetic field where 
energetic photons are emitted by electrons or positrons. In the supercritical 
magnetic field, the probability for a single photon decay into a pair can be 
approximated by a step function and thus, the free path corresponds to that 
for the pair production threshold. 

Since the outer gap is generally considered more efficient than the inner gap 
in conversion of the spin-down power to particle acceleration, we re-examine 
the possibility of neutrino production due to acceleration in the outer magne- 
tosphere. Ion acceleration in the outer vacuum gap, as a source of high energy 
neutrinos, was discussed long ago by Sato (1978) and recently by Bednarek 
& Protheroc (1997) and, as high energy cosmic ray sources, was considered 
by Bednarek & Protheroe (2002). In those cases, the detailed mechanism of 
ion injection and constraints by pair production on the acceleration were not 
discussed. In contrast to the inner gap scenario, acceleration of ions in the 
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outer gap, which is detached from the star's surface, requires a specific in- 
jection mechanism. Furthermore, the usual vacuum gap model (e.g. Cheng et 
al 1986) does not permit an external particle flux flowing through the gap. 
The effect of such external flux on the outer gap was considered recently by 
Hirotani & Shibata (1999), who showed that the gap location can be strongly 
modified due to injection of particles into the gap. Here, we consider the possi- 
bility that ion injection is due to a global current that circulate the pulsar- wind 
system (e.g. Shibata 1991; Hirotani & Shibata 1999) and that these ions are 
accelerated to ultra- high energies in the outer gap, leading to photomeson pro- 
cesses that produce high energy neutrinos. As for the polar gap scenario, one 
concentrates on the energetics and the particle flux, which can be determined 
respectively by the radiation-reaction limit and the relevant free path for pair 
production. The result should not depend much on the speciflc detail of the 
outer gap model. 

In Sec. 2, pair production processes in the MSM and constraints on the ac- 
celeration efficiency near the PC arc considered. Photomeson production on 
thermal radiation from the neutron star's surface is discussed in Sec. 3. The 
relevant neutrino flux originating from the inner magnetosphere is estimated 
and its detectability is discussed. In Sec. 4, proton acceleration in the outer 
gap is discussed. Photopion production in the outer magnetospheric region is 
considered in Sec. 5. The conclusions are summarized in Sec. 6. 



2 Acceleration of ions near the PC 



A neutron star forms with a hot surface on which free emission of charged 
particles occurs. We assume that for pulsars with fj,^^ -Q, < (where fi^ is the 
magnetic pole and Q is the angular velocity), free emission of ions (nuclei) from 
the PC occurs. These ions are accelerated along the fleld lines to ultra-high 
energies. Thermal photons from the star's hot surface can be Lorentz boost- 
ered to ultrarelativistic energy in the nucleus rest frame and produce pairs 
in the nucleus Coulomb field. Some electrons will be accelerated downward 
to initiate a cascade producing more pairs near the star, which provides seed 
positrons. These seed positrons are accelerated outward, initiating a further 
cascade above the PC. Whether the constraint on acceleration is determined 
by the proton free path or positron free path depends on nature of the polar 
gap. If the gap is nonstationary, the cascade by the seed positrons may play 
the more important role and if the acceleration proceeds in a steady, time- 
independent manner, the constraint is provided by protons. Both cases are 
considered here with particular emphasis on the former. 



4 



2. 1 Accelerating potential 



The maximum available potential across the PC can be written as 
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where Bq is the magnetic field on the PC, Bd = {2'kRq/ cPf/'^ ^ 0.46(1 ms/P)^/^ 
is the half-opening angle of the PC, Rq = 10^ cm is the star's radius. The spe- 
cific form of the potential depends on the detail of a particular model (e.g. 
Arons 1983; Harding & Mushmov 1998; Shibata et al. 1998). The basic as- 
sumption that is commonly used is that a nearly-parallel electric field develops 
due to the deviation of the charge number density of the outfiow from the GJ 
density ugj ~ B/2ZecP, where Ze is the charge of the particle. In the case of 
outflowing ions, since its composition is not well constrained by observations 
(Lai 2001), protons are considered. 

For acceleration sufficiently near the PC, the potential can be assumed to be 
a simple quadratic function of the distance from the PC (in units of Rq), 



where we consider the region s < 6d and e is the parameter that incorporates 
geometric effects such as the inclination angle between the rotation axis and 
the magnetic pole and the latitudes of the open field lines, the general relativ- 
ity effects such as the frame-dragging (e.g. Harding & Muslimov 1998, 2001; 
Hibschman & Arons 2001), and the detail how the low-altitude potential is 
matched to that at higher altitudes. For example, for a vacuum gap one has 
e ~ 1/0"^, which implies that the maximum potential is achieved at s ~ 
(Rudcrman & Sutherland 1975). If the potential is due to the space-charge 
effect, it increases with the distance more slowly than in the vacuum gap, 
with e ~ l/0d if the frame-dragging is dominant (e.g. Harding & Muslimov 
2001; Hibschman Sz Arons 2001) and e ~ 1 if the inchnation angle is nearly 
orthogonal (e.g. Harding & Muslimov 2001). Let (f)f = (p^es'j, where Sf is the 
characteristic distance where one pair is produced per primary particle. The 
efficiency can be estimated as 



provided that the energy loss is not dominant over the acceleration. It can be 
shown that only a very small fraction of the spin-down power can be extracted 
through particle acceleration in the magnetosphere of a young magnetar. 



= 0mes^, 
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2.2 The acceleration efficiency limited by pair production 



To estimate the acceleration efficiency limited by pair production, one derives 
the characteristic distance at which one pair is produced per primary protons 

or positrons (or electrons). First consider pair production by protons. Ultra- 
relativistic protons can produce pairs on a thermal radiation field through 
interaction of thermal photons in the proton's Coulomb field at a rate 

— ap^riphRo, (4) 



where ~ (3a/(TT/8vr) ln(4et/i7) ~ 5 x 10"^'' cm^ is the pair production 
cross section due to proton-photon collision (Chodorowski, Zdziarski, & Sikora 
1992), af ^ 1/137 is the fine constant, ar ~ 6.7 x lO^^^cm^ is the Thomson 
cross section, and Sth — 2.8© is the thermal photon energy in terms of the 
dimensionlcss surface temperature O = kBTs/mec"^ = 8.4x10"^ (Ts/5xlO'''K). 
We assume 2eth'y ^ 1- The photon number density at a radial distance r is 



nph ~ 0.36 



cks V r 

«.6.2xl0^3^^\V 7. j ^^_3^ 



r J V5 X lO^K 

where ass is the Stefan-Boltzmann constant. The proton free path can be 
estimated from 

(6) 



One estimates Sf ~ 3.2 x 10~^ for Ts = 5 x 10^ K, which increases slowly 
with an increasing pulse period P and rapidly with a decreasing temperature 
Tg. The proton free path is generally longer than that for an electron (or a 
positron) through single photon decay in a strong magnetic field. The latter 
is derived below. 

In deriving the electron (or positron) free path, it is worth noting two impor- 
tant features of pair production through a single photon decay in a supercrit- 
ical magnetic field B ^ Be- First, the probability for single photon decay into 
a pair can be approximated by the step function, i.e. one pair is produced 
when the photon energy reaches the threshold. Second, in the supercritical 
magnetic field, pairs are produced mostly in the ground state (Daugherty & 
Harding 1983; Weise & Melrose 2002), which tends to suppress further pair 
creation. The second feature implies that due to the absence of synchrotron 
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photons from the secondaries the pair cascade is less efficient than in moder- 
ately strong magnetic fields {B ^ Be). The first feature implies that estimate 
of the free path is reduced to calculation of the characteristic path for the pair 
creation threshold, given by (Harding, Baring & Gonthier 1997) 



> — — - — . for the || photon, (7) 

sm ' 

1 + (1 + 2£ Y^'^ 

£^> : — , for the ± photon, (8) 

where all relevant energies are dimensionless in units of mgC^, sb = B/B^ is the 
cyclotron energy, the _L and || modes correspond respectively to photons with 
polarization perpendicular to and in the plane defined by the wave vector k 
and the magnetic field B, and O^b is the propagation angle (with respect to B). 
Note that the threshold for the _L photon is higher than that for the || photon. 
Since in the superstrong magnetic field the ± photons can be efficiently split 
into the || photons, only the first threshold condition is relevant here. 

Two main radiation processes considered here that can produce pairs are res- 
onant inverse Compton scattering (RICS) (Herold 1979; Gonthier et al. 2000) 
and curvature radiation. RIGS occurs at ~ CthJ in the rest frame. The 
energy of the RIGS photon emitted by an electron (or positron) is ~ £^7 
(in the observer's frame). The || photon emitted at Sj will be absorbed and 
convert to a pair at s/, provided that 

Sf-Si> , (9) 



where = 4/3^^* > 4:/39d is the curvature radius (in Rq) of the field line with 
the magnetic colatitude 9^ on the PG. One has = A/Wd ~ 2.9(P/1 ms)^/^ 
for the last open field lines. The condition (9) can be derived directly from (7). 
For B ^ Be, the effective cross section for resonance scattering can be written 
as (Jeff = {3TraT/4:af)fB, where /s describes the relativistic effect that leads 
to the suppression of the cross section in the superstrong magnetic field. The 
photon production rate can be estimated from dNs/ds ^ Nm{'^m/l)'^, where 
= [9xQeB/{87^^1m)]i^ef[Ro/^l), K = h/ 171^0 is the Gompton wavelength, 
X — — ln[l — exp(— £5/97(1 — Pcos9m))], 9m is the maximum propagation 
angle of the incoming photon. For thermal radiation from the whole surface, 
one has cos^^^ = [(s + 1)^ — l]^/^/(s + 1). For Bq = 10^^ G, one may as- 
sume Jb ~ 0.01 (Gonthier et al. 2000). The photon production rate depends 
strongly on the distance s through 7 = 7f„es^, where 7^^ = e^^/mgC^ ~ 
1.3 X 10^^(1 ms/P)^(i?o/10^^ G). Therefore, one may obtain the following con- 
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dition for one pair-producing photon to be produced per primary: 



For Si <^ Sf, one may obtain Si ~ {Nm/'ie'^y^^ by neglecting 1/sj in (10). 
Prom (9), one derives 



10-3^1/3 / 5xlO^K y/^ /O^X V3 ._B^x -4/3 



j ■ (11) 



The gap length limited by RIGS is about s/ ~ 10^ cm/i?o for the above nomi- 
nated parameters. Since pairs are produced near threshold, in the ground state, 
further generations of pairs have to be produced by photons from RIGS or 
curvature radiation by the secondary pairs. The photons from the latter have 
energy too low to produce pairs. Since —E^n ~ 20^es//i?o ~ lO^^e Vcm~^ is 
much higher than that (2 x 10^°Vcm~^) required for ionizing positroniums 
(Usov & Melrose 1996), pairs produced inside the gap are mostly unbound. 
It is worth emphasizing that the characteristic length of a MSM polar gap is 
very small, much smaller than the PG radius. Thus, the quadratic potential 
(2) is a good approximation. Since a large backflow of particles may occur, 
the gap may be highly nonstationary with a characteristic frequency in the 
radio band, flc « c/sfRo ^3 x 10^s~\ 

Radiation from nonresonant inverse Gompton scattering may contribute to 
pair production. It can be shown that this process is less important than RIGS. 
In the Klein-Nishina regime, which is relevant here, the photon production 
rate is dNg/dt ~ —7/7, with 7 the energy loss rate. Since 7 ~ and —7 ~ 
(7rVTc/120A3)e2 ln(47e) ^ 2x W^{Ts/5 x 10^ K)^ ln(47e), one may estimate 
Si ~ (-7i?o/7mec). Prom (7) one finds Sf ~ aj^ms] ~ l/[ln(47msf©)]2. 

Pair production by curvature radiation can compete with RIGS in limiting 
the acceleration. The power of curvature radiation is 

P = 2^^^4 

-< curv 3 I ' K^'^J 



which gives the characteristic energy Ec ~ {\c/ Rc)'^'^ , where R^ = a^Ro- The 
production rate of curvature photons is Pcurv/ ^cf^eC^ ) which can be written 
into the approximate expression dNg/ds ~ (2re7m/3acAc)es^, where ~ 
e^/rUeC^ ~ 2.8 x 10~^^cm is the classical electron radius. Assuming at least 
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one photon in the energy Ec be emitted through curvature radiation, one may 
obtain Sf from Ng ~ (2re7r„/9acAc)e s^. The condition for A^^ = 1 yields 



The Lorentz factor of a positron at s/ is 7/ = 'jm^s'f ~ 5 x lO'^e^^^ for a 
milhsecond period, which is close to that limited by the energy loss due to 
curvature radiation, given by 7c ~ a^^/''7^/^ (3i?o/re)^'^^ ~ 8.2 x 10'''. The energy 
of the pair-producing photon is given by 



The photon must travel a further distance. As, before decaying into a pair. 
This distance can be determined from the threshold condition: Sc = 2ac/As. 
One estimates As ^ac/Sc ~ 10"^ ^ ^f- Therefore, pair production due 
to curvature radiation puts a much more severe constraint on the accelera- 
tion than RICS. The corresponding efficiency limited by pair production by 
curvature radiation can be written as 

^ , fj>^Y' ^ 3 , io-,,V3 (-1-)"" (^f. (15) 
' l2re7me/ UO^^GJ Urns J ^ ^ 



The efficiency limited by pair production by protons is given by 



V 



(cTp^nph-Ro)^ v5 X 10'' K 



-10-MfT74^ • (16) 



Other pair production processes such as pair production by photon-photon 
collision may also limit the acceleration near the pulsar. The pair production 
rate via 77 —> 6=*=, however, is significantly reduced in the superstrong magnetic 
field (Daugherty & Bussard 1980; Baring & Harding 1992). Since 77 oc e^^, 
the efficiency (15) is not very sensitive to the specific model of the potential 
and all the three cases discussed (cf. Eq. 2) lead to a very similar result except 
when the period increases substantially (cf. Sec. 2.3). 

In general, the efficiency determined by the proton free path is significantly 
higher (cf. Eq 6 and 16) than that determined by pair-production free paths 
of positrons (electrons). However, the latter is favoured as there seems no evi- 
dence from radio observations for two different classes of pulsars: one class with 
a proton-controlled gap and the other with an electron/positron-controlled 
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Fig. 1. Pulsar period as a function of time. The solid and dashed lines correspond 
to Bo = 10^5 G and 10^^ G, respectively. 

gap. In realistic situations, the polar gap acceleration is likely to be nonsta- 
tionary with protons playing only a minor role in the gap dynamics, and the 
gap energetics is limited by pair production by electrons/positrons. 



2.3 Evolution of the MSM polar gap 



Both and P arc time-dependent, which leads to time-dependence of the 
efficiency. In the case where pair production by electrons/positrons involves 
RICS or proton-photon coUision, the relevant free path depends on the surface 
temperature, which evolves with time. To derive the evolution of the MSM 
polar gap, we assume that the pulsar spin-down is due to the magnetic dipole 
radiation. Spin-down through gravitational waves may be important at the 
initial stage (e.g. Arons 2003) and such complication is not considered here. 
The pulsar period as a function of time is then given by 

P{t) = ^0 (l + -) , (17) 

where Pq = 1ms is the initial period, and = 'ic^ I B'IRq{Pq / 2nY ^ 2 x 10"^ s 
for Bq = 10^^ G and / = 10"^^ g cm^ is the typical spin-down time. The pulsar 
period P as a function of t is shown in Figure 1. The free path that defines 
the gap length is shown in Figure 2 for the cases of curvature radiation (solid 
line) and RICS (dashed line). The dash-dot lines represent the proton free 
path for Tg = 2 X 10^ K (upper) and 5 x 10^ K (lower). Since the spin-down 
time is much shorter than the cooling time of the neutron star, decreasing 
surface temperature does not modify much RICS. However, it does affect the 
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proton free path. Since the initial coohng processes are not well understood, to 
include this effect we model the star's cooling as a simple power-law function 
of time, i.e. Ts{t) ~ To(t/ls)~'', with 5 = 1/30 being used in in Figure 2 for 
protons and RICS. 

For MSMs, the effects of the RICS and curvature radiation on are compa- 
rable at the period 

This is in contrast with normal young pulsars in which RICS is generally the 
more important in constraining the energetics of the polar gap. Since in a 
supercritical magnetic field, RICS is suppressed and acceleration is extremely 
rapid (for P ~ 1ms), pair production due to curvature radiation is dominant 
over RICS until the period reaches the period (18). The efficiency r] as a, 
function of time is shown in Figure 3. For e ~ 1, this period is found to be 
Pc ~ 32 ms (cf. Figure 3) when t ~ 10®'^ s (about 46 day). The gap efficiency, 
initially constrained by pair production due to curvature radiation, increases 
as the MSM spins down. When the RICS becomes dominant at Pc, r] decreases 
gradually because effective RICS favours moderate acceleration. Generally, the 
efficiency decreases to a minimum and then the effect of the surface cooling 
takes over, which renders an increase in rj due to a decreasing Tg. 



3 Photomeson production 

The outward-propagating relativistic protons can produce neutrinos on the 
thermal radiation from the star's surface through the photomeson process. The 
dominant channel for neutrino production is photomeson production at the A 
resonance (Stecker et al. 1991; Waxman & Bahcall 1997), corresponding to the 
process 7p — > A — > mr'^. The pion decay, tt"*" — > //"^ -|- i/^ — > e"^ -|- + + f^, 
produces neutrinos. For the photon energy e'^ in the rest frame of the proton, 
the cross section is peaked at e'^ ^ Eq = 0.35/ip with the width A^o ~ 0.12/ip 
(Stecker et al. 1991). The Lorcntz factor 7^ = 1/(1 — Pl)^^"^ of the proton at 
which the cross section for pion production peaks is estimated to be 

' {1- (3, COS 9)eth V Jl-(3,cos9' ^ ' 

where 9 is the propagation angle of the incoming photon, given by cos^ = 
[(1 -\- sY — l]^/^/(l + s). Since one has cos6' pa (2s)^/^ for s -C 1 near the 
surface, the factor 1/(1-/3* cos 9) can be discarded. Although the proton free 
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Fig. 2. The free path for pair creation. The dashed hnes is for RICS, which from top 
to bottom correspond to e = 1/^^, l/^rf and 1. The three soHd Uncs is for cTirvature 
radiation, corresponding to e = 1/^^ (lower), l/6d (middle) and 1 (upper). One 
assumes ~ 5 x 10^ K and Bq ~ 10^^ G. For comparison, the proton's free path is 
shown as the dash-dot lines for = 2 x 10^ K (upper) and 5 x 10^ K (lower). For 
e = 1, RICS takes over the effect of curvature emission at about t ^ day when the 
period increases to about P ~ 32 ms. For other two cases, RICS becomes dominant 
at t > 10^ s (not shown in the figure). Since 6d oc and 6d ~ 5.5 x 10~^ at 

t ~ 10^ s, one always has s / ^ in the parameter regime considered here and the 
quadratic potential (2) should be a good approximation. 

path gives higher 77 than electrons/positrons (cf. Figure 3), which may be 
applicable for a steady gap, wc emphasize the latter case which is applicable 
when acceleration is nonstationary (cf. Sec. 2). Assuming the characteristic 
acceleration length is S/, which is limited by pair production by seed positrons, 
one obtains the maximum Lorentz factor of the accelerated proton, 

e0m^ Iran .^„x 

7/ = 2 = • (20) 



Photomeson processes at the A resonance requires that the Lorentz factor 
of the accelerated protons satisfies the resonance condition 7/ ~ 7*, which 
occurs near the peak of the thermal photon distribution. One can show from 
(15) that for all three cases considered (e = 1/9^, e — l/9d, and e = 1), the 
threshold 7 > 7*(1 — Aeo/^o) is satisfied. The pion energy is e,, ~ 0.2ep, where 
£p = 7*A*p is the proton energy (in mgC^). Assuming that the energy is evenly 
distributed among the four products and that pions decay before significant 
cooling or acceleration occurs, one may estimate the energy of neutrinos as 
£1/ = £7r/4 ~ 0.05£p lO^Hp, about 1 TeV. 

In a strong thermal radiation field, energy loss of charged pions due to in- 
verse Compton scattering can be important. The decay time is — 2.6 x 
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Fig. 3. The evolution of the efficiency for r^(0) = 2 x lO'^ K (dashed) and 5 x 10^ K 
(soUd). In each case, Hnes from top to bottom correspond to e = 1/0^, 1/^^ and 1. 
The break corresponds to that pair production due to curvature radiation switches 
over to that due to RICS for e = 1. For the other two cases, the breaks occur at 
a later time (not shown in this figure). The dash-dot-dashed line corresponds to 
the efficiency determined by the proton free path (for Ts{0) = 5 x 10^ K), which 
is significantly higher than that by positrons/electrons. As in Figure 2, we assume 
Bo = 10^^ G. 

10~*77rS ^ 5 X 10~^s for = 3.7 x 10*^ in the observer's frame. The energy 
loss rate of pions is ~ caT-KSthnphl'iijne/mT,), where the cross section is 
cttt = {nT-e/m-^yaT ~ 1-3 x 10~^o't, 171.^ ~ 274me is the pion mass (about 140 
MeV; and the neutral pion mass is about 135 MeV). The coohng time can be 
obtained from 



tic ^ 



10" 



X 10^ 



r 



-2 / 3.7xl0^ ' 



(21) 



Since the threshold for pion production is reached before the proton reachs s/, 
charged pions produced by the proton can be accelerated and this acceleration 
time is much shorter than either the cooling time and decay time. The pion has 
been accelerated to very high energy before it decays. Therefore, the typical 
pion energy is 8^, = 0.27*/ip + (7/ - 7*)Aip = (7/ - 0.87*)/ip ~ 10'', which would 
give rise to the neutrino energy e^, ~ £^/4 ~ 2.5 x 10^ or 1.3 TeV. However, 
in the region above Sf, the pion suffers energy loss due to IC and the the 
neutrino energy can be much lower than the above estimate. The condition 
= ti^ leads to 7^ = 1.2 x 10^(5 x WK/T^f at r ~ i?o and 7^ ^ 5.6 x 10=^ 
at r ~ Rlc, corresponding to Sy = Et^/A pa 8.2 x 10^ (or 40 GeV) and 3.8 x 10^ 
(or 2 X 10^ GeV). 
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The optical depth for photomeson production can be approximated by r^- ~ 
n-pi^a^RoAs, where 0",^ ~ 2 x 10^^^ cm^ is the cross section at the A resonance 
peak, As is the characteristic length for pion production, and riph is the photon 
number density given by (5). Assuming = 1 one obtains As ^ l/up^^aT^Ro ~ 
0.12 at r Rlc — 4.7 x 10^ cm and one would have a much smaller As 
than this if r ~ Rq. The neutrino cmissivity is A„ ~ [e.,^ /4)a.jrnphnGjC with 
Ugj ~ 3.5 X 10^^(-Bo/10^^ G)(P/1 ms)^^ cm^'\ from which one may derive the 
neutrino luminosity L^, ~ 4:7T{9dRo)^RoA^. The neutrino number flux density 
can be approximated by 



AQoD^e^ ■ AQoVlmsy \10^^G 



-2 



x( Il—)'(^L_] cm-^s-\ (22) 

\5xWkJ llOMpci ^ ^ 



where AQq is the beaming angular width, D is the distance to the source, and 
one assumes AsPo'^vr^ph = 1- Since neutrinos are generally detected through 
production of muons, the event rate can be calculated from N^, ^ F^P^^, 
where P„i^ = 10"'^(£i,/2 x 10^), with = 2 x 10^ corresponding to 1 GeV, is 
the probability of the muon event (Halzen & Hooper 2002). Then, one finds 



AT KA ( P \f 



AQo Vims/ VlO^^G/ V5 X lO^K 



^UoMpcJ V 100 GeV 



km^Vr"^ (23) 



The event rate in the early spin-down phase is about A^i^ = 54km~^yr~^ 
with the neutrino energy about 100 GeV. Such rate may be too low to be 
detectable with the current AMANDA and requires a km scale detector with 
a lower energy threshold of about 100 GeV. 



The photomeson process also produces relativistic neutrons, which move bal- 
listically through the magnetosphere and may interact with target nuclei in the 
stellar envelope to produce neutrinos and gamma-rays (Protheroe, Bednarek 
& Luo 1998). The neutron decay time is i„ ~ (920s)7„, where 7„ is the 
Lorentz factor of neutrons. For 7n ~ 7/ ~ 5 x 10^, the decay distance is 
about t„c ~ 1.3 X 10^^ cm, much larger than the size of the shell ~ f^t, where 
Vg ~ 3000 km s^^ is the expansion speed. If the shell is completely fragmented 
(Arons 2003), neutrons would propagate freely through it and decay outside 
the remnant. 
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4 Proton acceleration in the outer magnetosphere 

Protons may be accelerated to much higher energy near the LC where screen- 
ing of the electric field occurs in the perpendicular (to the field line) direction 
and the maximum energy is limited only by radiation loss. To accelerate pro- 
tons in the outer magnetosphere, a specific mechanism for injection of protons 
into the gap is needed and will be discussed in Sec. 4.1. 



4.1 MSM outer gap 

In close analogy to the outer gap in normal young pulsars, we assume an 
oblique rotator and that an outer gap develops in a region from the null surface 
along the open field lines to the LC, with an upper boundary surface defined 
by a pair production front and a lower boundary by the last open field lines. 
As in the case of the polar gap, such gap can be nonstationary due to the very 
nature of time-dependent pair cascades. The gap energetics is thus determined 
by the average thickness (between the two boundary surfaces) Wg of the gap, 
where Wg is determined as in the PC scenario by the pair production free path 
of an electron or a positron (see discussion below Eq 24). Since the available 
potential Wg(f)m well exceeds that for radiation-reaction limit of protons due 
to curvature radiation, the oscillatory nature of the gap can cause modulation 
of the particle fiux and it does not constrain the proton's maximum energy 
(it is limited by radiation reaction). One may estimate the maximum proton 
energy and relevant particle flux by following similar procedure to that used 
in the steady outer gap model (e.g. Cheng, Ho, Ruderman 1986; Romani 
1996; Hirotani, Harding & Shibata 2003) and the result should not depend on 
particular aspects of the model. 

It is assumed that particles are accelerated along the fleld lines to ultra-high 
energies and emit high energy photons through curvature radiation or inverse 
Compton scattering and these photons convert into pairs through single pho- 
ton decay in the strong magnetic field. Note that in normal young pulsars the 
only effective channel for pair creation near the LC is through photon-photon 
collisions. Since the field lines curve away from the magnetic axis, the pairs 
are produced on upper, neighbouring field lines and thus one may estimate 
Wg in a way similar to that for the polar gap, i.e. by estimating the pair pro- 
duction free path of an electron (or a positron). Because the magnetic field 
near the LC is weaker than near the PC, the pair production free path should 
be calculated from the opacity of photon absorption (due to pair production), 
characterized by a parameter x = (l/2)£7eB sin^fc^, where e.y is the energy of 
pair producing photons and 9kB is the photon propagation angle as defined 
in (7) and (8). Apart from the threshold conditions (7) and (8), the foUow- 
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ing condition is needed for producing one pair (corresponding to the opacity 
T ~ 1) (e.g. Arons 1983): 

X > A f^) . (24) 



15 VlnA 



where the near-threshold effect on the opacity (Daugherty & Harding 1983) 
is ignored and In A is a logarithmic parameter defined through A(lnA)^ — 
1.6{af/e'^6B){Rc/K), which is not sensitive to the pulsar parameters, with a 
value between 1 and 30. For — 5 x 10^, one has InA ~ 2. For a photon of a 
given energy e^, the typical angle at which the photon converts into a pair is 

si„«.. . 1.2 X 10-r {^f {^f . (25) 



where ^ = v/Rlc- The characteristic thickness (in Rlc) of the gap can be 
estimated from Wg ~ sin^fc^. The time dependence of the LC radius and the 
magnetic field at the LC are shown in Figures 4 and 5, respectively. Single 
photon decay should be the dominant pair production process up to t ~ 10^'® s 
when the magnetic field at the LC drops below 10^ G. When P is sufficiently 
long, magnetic pair production is no longer important and the only channel 
for pair production is through photon-photon collisions and like normal young 
pulsars the MSM becomes more and more efficient with increasing P to a 
maximum where the outer gap is completely open. 

One problem with acceleration of protons in the outer magnetospheric region 
is the proton injection. An outer gap model including injection of particles 
was recently discussed in detail by Hirotani & Shibata (1999), also Hirotani, 
Harding & Shibata (2003). (In their discussions, protons are not included in 
the model.) It has been recognized that current circulation in the pulsar- wind 
system may play an essential role in transforming rotational energy to particle 
kinetic energy (e.g. Shibata 1991). Such global current should exist even in 
the case of nonstationary acceleration. Thus, one may assume that the ion 
injection is due to the global current system. The sign of the accelerating 
electric field is determined by the sign of B^. For an outer gap located outside 
the null surface {B^ — 0), the accelerating field is positively directed when 
Bz < 0, corresponding to pulsars with /x„ ■ fi > 0. Therefore, outflowing 
primary electrons from the PC and downflowing positrons (possible protons) 
form a return current directed toward the PC, and a flow of protons plus 
positrons through the outer gap forms part of the outward-directed current. In 
the steady model, inclusion of an external flux shifts the location of the zero G J 
density, the 'null surface', and hence the gap location. For a flux injected from 
the null surface, the zero GJ density is shifted towards to the LC, which shifts 
the gap to the LC (Hirotani, Harding & Shibata 2003). Although such features 
are derived based on the steady assumption, they may well be applicable to 
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the nonstationary case if acceleration confined to the region very close to the 
last open field lines where the influence of nonstationary pair cascades is the 
least. 



4-2 Maximum proton energy 



The outer gap potential, (pm^g, limited by pair production by positrons/electrons, 
is much larger than that of the polar gap, (pms}- For example, at P = 1ms, 
the outer gap has lO~^0m {wg ~ 10^^ at F ~ 1 ms), as compared to the polar 
gap 10~^(f)rn (s/ ~ 10^'^, cf. Eq. 13). For the relevant parameter regime one 
is interested in, pair cascades do not constrain the acceleration length since 
pairs are produced on field lines with a smaller magnetic colatitude. So, the 
maximum proton energy is limited only by radiation reaction such as curva- 
ture radiation. Since the specific form of the potential is model dependent, 
here we write the parallel accelerating electric field in the form 

^11 = (26) 



where /„ < 1 is the fraction of the full (vacuum) potential across the gap, 
tns/Rlc < ^ < 1, tns ~ -Rlc(2 cot CKj/S)^ (in the Q,-B plane) is the radial 
distance to the null surface, is the MSM inclination angle (relative to the 
spin axis), and 0™, is given by (1). The power index is in the range of —3 < 
a < 0. In Romani (1996)'s model, a = — 1 is assumed. For an oblique rotator. 
Tats is a substantial fraction of the LC radius, one may use an approximation 
^~1. 

The maximum proton Lorcntz factor can be derived by equating the energy 
loss rate due to curvature radiation to the acceleration rate. Note that the 
power of curvature emission for protons is the same as that for electrons, 
given by (12) with 7 now being interpreted as the proton Lorentz factor. The 
radiation-reaction limited Lorentz factor is derived as 



/ ?>-nn.faRLc \ ^'^ ( Wg_\ (2„+5)/8 

= r.6xioVi/^(^)"'\^)'^(^)'^^-v«. (27) 

For fa = 10-2, one has 7c ~ 2 x 10^{wg/^^/'^y/\P/lmsy/^. The radius of 
the open field line region can be written as R± ~ ^^^^Rlc- The condition 
Wg < R±/Rlc implies that Wg/^^/'^ < 1. As one is interested only in the outer 
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Fig. 4. Evolution of the light cylinder radius. The solid and dashed lines correspond 
respectively to Bq = G and 10^^ G. 

log(BLc/G) 
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Fig. 5. Evolution of the magnetic field Blc at the light cylinder radius. The solid 
and dashed lines correspond respectively to Bq = 10^^ G and 10^^ G. 

gap near the LC, ^ ~ 1, the radiation- reaction limited Lorentz factor is not 
sensitive to the specific choice of a value. 



5 Photomeson production neeir the LC 

In the outer magnetospheric region, both thermal radiation from the star's 
surface and nonthermal radiation from cascades due to particle acceleration 
in the outer gap can provide target photons for photomeson processes. Accel- 
erated protons may also interact with soft photons from the diffuse thermal 
radiation, e.g. from the hot ejecta or the interaction region between the MSM 
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wind and the remnant shell. Here, we consider the first and third possibilities 
only. 



5.1 Photomeson threshold 



The thermal radiation from the star's surface subtends a much smaller solid 
angle near the LC than near the PC and thus the protons satisfy the pho- 
tomeson threshold at much higher energies. Assuming that the acceleration 
occurs near the last open field lines, the characteristic propagation angle of 
thermal photons originating from the star's surface can be derived from 

cosOm^l-^, (28) 



where <^ 1. Since this condition is generally satisfied, the geometric effect 

of the thermal radiation is important, which leads to the following expression 
for the photomeson production threshold for the photon energy, 

Sph ~ 0.35^. (29) 



Since ^'^/O^ ^ Ij protons satisfy the threshold at much higher energies than 
near the PC, producing pions with extremely high energy. Eq. (29) implies 
that Eph ~ P^, higher proton energy is required as the MSM spins down to a 
longer period. At P ~ 50 ms, one has 7 > 10^. 

Charged pions are subject to the energy loss due to IC. Since the Lorentz 
factor of pions is 7^ > 0.2'J^,{mp/mn), the photon energy in the pion rest frame 
ise' ~ 2.807^,^72^^ > 1.3xlO^(mp/m^)(5xlO^K/r,)-i > 1. In the following 
discussion, the surface temperature is assumed to be time independent. 
Hence, IC is in the KN regime with a characteristic time: 



~7r2(7Tce2^^1n(47^e) 

We assume that the thermal radiation field originating from the neutron star's 
surface has spherical symmetry with the photon number density decreasing 
radially as The LC radius increases as the MSM spins down (cf. Fig- 

ure 4). Therefore, the photon number density near the LC reduces as the 
magnetosphere expands. One may find the characteristic period Pi at which 
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the energy loss time of charged pions near the LC is equal to the pion decay 
time. This period is given by 

Pi ~ (11.8ms)(T,/5 X 10^ K) M^lfn^)V^\ ^3^^, 



Charged pions lose much of their energy through IC until the MSM spins down 
to P ~ Pi ~ 21 ms when the photon number density is sufficiently low and IC 
can no longer constrain the neutrino energy. The neutrino energy can reach 
- 0.057c/Xp = 5 X 10^, about 10^ TeV. 

If the region surrounding the MSM is filled with diffuse thermal radiation, 
extremely energetic protons produce pions on these soft photons. There are 
two possible sources for the diffuse thermal radiation: the hot ejecta in the 
early phase of the pulsar spin-down (e.g. Beall & Bednarek 2002), and the 
interaction region of the inner surface of the remnant where the MSM wind 
may deposit its energy (e.g. Rees &; Mcszaros 2000). The threshold condition 
for photopion production is 7 ~ 3.2 x 10^(5 x lO^K/T^j), where is the 
temperature of the thermal radiation. Here we assume that the radiation field 
is isotropic. This condition is well satisfied by (27). For such a low T^, pion 
cooling is not important in constraining the neutrino energy. 



5.2 Neutrino flux 



The neutrino fiux can be estimated from the proton luminosity, given by 



;iO-e.gs-)(^)"f^l^. (32) 



imsy Vioisc; ^3/2 ' 

where the proton fiux in the gap is assumed to be a fraction /p < 1 of the 
GJ fiux ugjc. The neutrino number fiux density can be estimated from F^, ~ 
L^/47rAQo^^£z. ~ (4.3 x lQ-^'^c^-'^s-^){P/l^s)-^{fpWg/^^''^) for AQq = 
47r, Pq = 10^^ G and P = 10 Mpc. The atmospheric background neutrinos 
have the fiux spectrum ~ 10~'^(e,y/2 x 10^)~^'^ cm~^ s^^ sr^^. The planned 
IceCube has an angular resolution of about 1°. The background number fiux 
density from a smaU patch (tt/ISO)^ of the sky is P^ ~ 0^(7r/18O)^ ~ 9.6 x 
lO-^'^ cm-2 at £^ ~ 2 X lO'^ or about lO^TeV. The predicted fiux is well 
above this background fiux for /p > 10~^. 

One estimates the maximum neutrino event rate as N^, ~ F^P^^ with P,^^ ~ 
2 X 10~^(£,,/2 X 10^) for neutrinos above TeV energies (Halzen & Hooper 2002), 
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leading to 

,8x loV.(^)" (^) (j^f k,n-yr- (33) 

where one assumes Bq = 10^^ G. For example, for P ~ 48 ms, for which one has 
Wg/^^/'^ ~ 0.13, one has an event rate of about Ny ~ (1.1 km~^ yr~^)(/p/Ar2o), 
which is potentially detectable by the IceCube provided that the proton flux 
in the gap is a substantial fraction of the GJ flux and that the emission is 
moderately beamed, say Afio ~ 1 and the source is relatively nearby. Since 
the specific value of /p < 1 is not well constrained, one cannot exclude the 
possibihty that the current that forms a closed circuit is predominantly due to 
positrons. In this case, the neutrino flux produced from the outer gap would 
be strongly hmited by the factor fp <^ 1. 



6 Discussion and conclusions 

Neutrino production in the magnetosphere of a fast rotating magnetar is con- 
sidered. It is shown that MSM can be a strong source of high energy neutrinos 
in its early spin-down phase. We consider proton acceleration in both the inner 
magnetosphere near the PC and the outer magnetosphere near the LC. When 
protons are accelerated near the PC, their maximum energy is constrained by 
pair production due to accelerated seed positrons/electrons or primary pro- 
tons. In supercritical magnetic fields, photon splitting may be important in 
limiting pair production (Baring & Harding 2001). However, since this pro- 
cess can only split the ± polarized photons (Baring & Harding 2001; Usov 
2002), it may not prohibit pair production by the || photons. Pair annihila- 
tion is greatly enhanced in a supercritical magnetic field leading to reduction 
in the secondary pair density, but such process does not strongly affect the 
electric field screening. Due to the supercritical magnetic field, the efficiency 
of neutrino production can be estimated from the pair production free paths 
of positrons/electrons. The gap determined by pair production by protons has 
a relatively higher efficiency (cf. Figure 3). Although the proton-controlled 
steady gap cannot be excluded, the assumption that electrons/positrons con- 
trol the dynamics of the polar gap appears to be the more consistent with the 
lack of any observational evidence for two classes of radio pulsars. If protons 
play a role in determining the polar gap, one would have one distinct class 
of pulsars with their gap controlled by electrons/positrons and the other by 
protons. For a typical radio pulsar with Ts ~ 10^ K, the proton-controlled gap 
is very inefficient in producing pair plasmas that needed for production of the 
observed coherent radio emission. So far, there is no obvious observational 
evidence in pulsar radio emission for this distinction. In the case of the gap 
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limited by pair production by electrons/positrons, photomcson production is 
possible provided that the surface temperature is T, ~ 5 x 10'' K. Because of 
the limit on the proton energy by pair production, only a tiny fraction (about 
7] ~ 10~^ initially) of the spin-down power goes into protons and efficiency 
increases as the MSM spins down while the proton flux decreases. As pho- 
tomeson production occurs near the PC, charged pions are subject to energy 
loss to inverse Compton scattering, which further limits the neutrino energy 
to about 100 GeV, below the threshold for IceCube. 

Acceleration in the outer gap can be more efficient than in the polar gap as 
the maximum proton energy, limited by curvature emission, can reach about 
lO^mpC^. It is suggested here that proton injection in the outer gap is part 
of a closed global current that flows through the gap in the form of proton 
flux. In the application to MSMs with fjt^ • 17 > 0, the outflowing particles in 
the PC region are mainly electrons that, together with possible downfiowing 
positrons or protons, provide a return current into the PC. The outer gap is 
limited by pair production due to magnetic single photon decay. The efficiency 
is initially small due to efficient pair production limiting the gap thickness, 
which in turn limits the total flux of protons that can pass through and ac- 
celerated in the gap. The efficiency increases with increasing period. Since the 
thermal radiation from the surface subtends a much smaller solid angle, pro- 
tons satisfy the photomeson threshold at much higher energies so that pions 
are produced at high energies ~ 10 TeV. As the star spins down rapidly, the 
magnetosphere expands and the photon number density in the outer magne- 
tosphere decreases. When P ~ Pi, charged pions decay into neutrinos before 
they lose their energy and TeV neutrinos emerge. TeV neutrinos can also be 
produced as a result of interactions of ultrarelativistic protons originating from 
the magnetosphere with soft photons from the diffuse thermal radiation due 
to the hot ejecta, which may exist in the early phase of a new born MSM, 
or due to the heating of the interface region between the relativistic MSM 
wind and the remnant shell. The neutrino number flux density is estimated 
and it is shown that the corresponding muon event rate may be detectable 
with the planned IceCube provided that the source is relatively nearby or 
the emission is moderately beamed. The estimate is based on the assumption 
that the outer gap is time-independent, which may not be realistic. Inclusion 
of time-dependence of pair production and hence the nonstationary gap may 
lead to a signiflcant modiflcation to the particle injection and acceleration. 
A nonstationary gap does not have a well defined pair production front and 
hence a much larger effective cross section area, allowing more protons to be 
accelerated leading to an increase in the neutrino luminosity. In the above 
discussion, one ignores the possibility that a positron flux forms a major part 
of the current that closes the current circuit (e.g. Shibata 1991), with a proton 
flux limited to a small fraction of the GJ flux. If this is the case, the neutrino 
flux is severely constrained by the fraction factor fp <^ 1. 
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So far, the polar gap and outer gap scenarios have been treated separately. 
In realistic situations, pair cascades near the polar cap may affect the outer 
gap and vice versa. To determine self- consistently such complication requires 
a quantitative model that hnks the two acceleration regions and such model 
is currently not available. Apart from the polar gap and outer gap scenarios, 
protons may be accelerated to ultra-high energy in a slot gap (Arons 1983; 
Muslimov & Harding 2004). In such model, protons can be assumed to be 
primary particles extracted from the PC and can be accelerated to energies 
well exceeding the pion production threshold (e.g. Protheroe, Bednarek, & 
Luo 1998). Since a slot gap may accelerate protons to the maximum energy 
limited by radiation-reaction and the available particle flux (through the gap) 
is constrained by pair production, the result (the neutrino energy and flux) 
discussed in Sec. 5 should be qualitatively valid for this case as well. However, 
a quantitative prediction of the neutrino flux from the slot gap requires further 
work. 
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